Radio Properties of the Supernova Remnant N157B 
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ABSTRACT 

A new investigation of the supernova remnant (SNR) N157B was carried out with the Aus- 
tralia Telescope Compact Array. Radio continuum images of the entire 30 Doradus region have 
been made at 3.5 and 6 cm wavelength with a resolution of ~ 2". These data allow a high res- 
olution study of the spectral index distribution and polarization properties of both N157B and 
the nearby 30 Doradus nebula (the latter will be reported in a subsequent paper). N157B is an 
extended Crab-type SNR which may be beginning the transition to a composite remnant. There 
is little apparent fine structure and the brightest radio region is several parsecs from the probable 
position of the X-ray pulsar. The SNR has a radio spectral index of —0.19 and is significantly 
polarized at 3.5 cm but not at longer wavelengths. 



Subject headings: radio continuum: interstellar - 
individual (LMC) 

1. Introduction 

N157B (Henize 1956) (also 30 Dor B or SNR 
0538-691) is a supernova remnant located near the 
edge of the 30 Doradus nebula (Bode 1801) in the 
Large Magellanic Cloud (LMC). The 30 Doradus 
complex is the nearest extragalactic giant H II re- 
gion and an active star-forming region (Chu et 
al. 1992). The nebula and the SNR are in the 
same interstellar region of the LMC as SN1987A 
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and other SNRs, among which is SNR 0540-693, 
a composite remnant with a pulsar and a plerion 
component (Manchester et al. 1993). The SNR na- 
ture of N157B was first determined by its excess of 
radio over optical emission (Le Marne 1968), and 
confirmed by the detection of [S II] lines in the 
optical spectrum (Danzinger et al. 1981). 

The Crab nebula represents the class of pleronic 
supernova remnants; it contains a pulsar as a cen- 
tral energy source and a well-organised magnetic 
field indicated by strong linear polarization at ra- 
dio wavelengths and relatively flat radio spectrum 
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(Weiler 1983). A number of detected SNRs have 
radio characteristics similar to the Crab nebula, 
yet the pulsar stimulation of the synchrotron neb- 
ula is confirmed for just a few of them (Seward 
1983, 1989; Weiler 1988). N157B is classified as 
a Crab-type remnant because of its centre-filled 
morphology, flat radio spectrum (Mills et al. 1978; 
Milne ct al. 1980), and strong soft X-ray emission 
(Long & Hclfand 1979). 

Recently a 16 ms pulsar detected only in X- 
rays (Marshall et al. 1998; Gotthclf & Wang 1996; 
Wang & Gotthelf 1998) has been associated with 
N157B. Previous observations of N157B at radio 
wavelength have not detected any polarized emis- 
sion (Dickel et al. 1994). However, our current in- 
vestigation at a shorter wavelength shows a signif- 
icant level of polarization. All of the data together 
confirm the classification of N157B as a Crab-type 
SNR. 

The presence of N157B in a relatively dense and 
complex interstellar environment together with 
the younger SN1987A and the presumably older 
SNR 0540-693, offers a unique opportunity to 
study the evolution of Crab-type SNRs in such 
areas. The best Galactic examples, the Crab and 
3C58, both appear to be in exceptionally low- 
density environments (Fesen 1997). 

2. Observations and data reduction 
2.1. General procedures 

The first obervations at 6 cm were carried out 
with the Australia Telescope Compact Array 2 
(ATCA) in the period January - July 1993. Addi- 
tional observations at 3.5 and 6 cm were obtained 
during May 1998 - June 1999. The six antennas 
were used in several different array configurations 
giving baselines with the range from 31 m to 6 km 
to provide good sampling of the observed region. 
Details of the ATCA can be found in Frater et 
al. (1992). The relevant observing parameters are 
given in Table 1. 

The primary calibrator used in all observations 
was PKS B1934-638 with assumed flux densities 
of ^6 and 2.78 Jy at 6 and 3.5 cm respectively. 
PKS B0454-810 was used as secondary calibration 

2 The Australia Telescope is funded by the Commonwealth 
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by CSIRO 



source. 

The observations in 1998 - 1999 were done in 
mosaic mode with 20 pointing centres to cover the 
whole 30 Doradus region. The observing pattern 
was such that the pointing nearest to the SNR 
(field no. 15 - see Table 1) was observed after ev- 
ery five fields to improve signal-to-noise on this 
object. The data were taken simultaneously in 
32 x 4 MHz frequency channels of which the 26 
central channels are used in the imaging proce- 
dure. They have been averaged down to 13 x 8- 
MHz channels by Hanning smoothing and then to 
a single channel during imaging. The data were 
gridded and deconvolved according to the proce- 
dure for mosaic data sets (Sault & Killeen 1997). 
Separate images of field no. 15 alone were also pro- 
duced. 

The 3.5-cm images of the SNR were made by 
combining 7 fields (including no. 15) which had the 
pointing centre closest to the SNR. The deconvo- 
lution of all 7 fields was done together, which im- 
proves the deconvolution process but introduces 
some errors in the model of the primary beam. 
The edges of the mosaiced region have low sen- 
sitivity so the noise level across the mosaic im- 
age will vary. In the imaging process the pri- 
mary beam attenuation will not be completely cor- 
rected, because that correction is constrained by 
the noise in the individual pointings, which can 
differ (Sault & Killeen 1997). In this case, the pro- 
cedure limits the dynamic range of the final image 
to 100:1. We also made an image of field no. 15 
alone. This was cleaned and corrected for the pat- 
tern of the primary beam. The peak brightness in 
both images is very close (~7mJy/beam), but the 
single field image has lower signal-to-noise ratio 
over the entire field. 

The image at 6 cm was constructed by combin- 
ing four data sets at 4.8 GHz and pointing no. 15 
from 4.4 GHz data (see Table 1). Thus, in addi- 
tion to the mosaicing procedure explained above, 
the multi-frequency synthesis technique (Sault & 
Killeen 1997) was also used to produce the final 
image of the region. As in the previous case, an 
image of field no. 15 alone was made as well. 

2.2. Polarization 

The interferometer produces an instrumental 
polarization which increases outward from the cen- 
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Table 1 
Observational parameters. 



Parameter 


4.4 GHz 


4.8 GHz 


8.6 GHz 


Date 


1998 Aug 27 


1993 Jan 16, 19, 29 


1998 May 1, 31, Aug 27 




1998 Sep 26, Oct 2 


1993 Jul 23 


1998 Sep 26, Oct 222 




1999 Jun 24 




1999 Jun 24 


Configuations 


6C, A, 1.5D, 375 


750B, C, A, D 


750A, E, 6C, A, 1.5D, 375 


Field Center (J2000) 


05' l 37 m 50 sa 


05 /l 38 m 45 s 


05 ft 37 m 50 ;sa 




-69°10'00" a 


-69°06'00" 


-69°10'00" a 


Total bandwith (MHz) 


128 


128 


128 


No. of frequency channels .. 


32 


32 


32 


FWHM of primary beam ... 


10' 


10' 


5' 



a field no. 15 of a larger mosaic 



ter of the primary beam pattern of the individual 
telescopes. In fields with the pointing centre away 
from the SNR for which an unpolarized source (30 
Doradus) is at the half-power point, we found frac- 
tional polarizations up to ~10% with electric field 
vectors radially oriented with respect to the di- 
rection of the beam centre. This value of frac- 
tional polarization in the surrounding fields con- 
firms that off-centre instrumental polarization is 
indeed strong. 

Therefore, for the polarimctry we used only 
field no. 15 for which N157B is in the center of 
the beam and the instrumental polarization is low 
(0.3%). The images of the Q and U Stokes param- 
eters were made in a similar way as the / image 
and combined to produce the image of polarized 
intensity. To create such an image we used the 
MIPJAD task IMPOL which statistically corrects 
the polarized intensity for Riccan bias (Sault & 
Killccn 1997). Noise in both the Q and U polar- 
ization images was below the 5-cr noise level of 
the total intensity image and CLEANing was not 
necessary. 

The polarized intensity images of the SNR at 
6 cm were made using only field no. 15 with the 
same procedure as for the single-field image ex- 
plained above. The 1993 positions were not di- 
rectly toward N157B and so were not included in 
the polarization analysis because of the off-axis 
effects disscussed above. The fractional polariza- 



tion in the direction of N157B was found to be 
~3%, near the instrumental value, thus confirm- 
ing only that N157B is not highly polarized at this 
wavelength. 

3. Results 

3.1. 3.5 cm image 

The high resolution total intensity image at 
3.5 cm of the region around N157B was made using 
all six available configurations; it has a Gaussian 
beam of HPBW ~2" (see Fig.l). It shows rea- 
sonably uniform brightness distribution and the 
center-filled morphology characteristic of Crab- 
type SNRs. We note that the bright peak is re- 
solved. It lies about 12" (4pc at the 50kpc dis- 
tance to the LMC) to the northwest of the X- 
ray peak, which is presumably the location of the 
pulsar (Wang & Gotthclf 1998). Also, the radio 
peak is not centered on the full SNR which ex- 
tends well to the south with dimensions of about 
120" x 70" or 30pcxl8pc. Slices along the major 
axis of the remnant at a position angle of about 
150° show continuous emission throughout the to- 
tal extent with a slight decrease near a declination 
of — 69°10'50" and then a rise again to the south. 
Perhaps part of an outer shell may be forming 
around the plerion. The size given includes this 
outer emission. 

This is the most sensitive high-resolution radio 
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image of this object now available, yet no direct 
evidence of the pulsar can be seen. This is not sur- 
prising, however, as 3.5 cm is a very short wave- 
length for pulsar emission which usually peaks at 
meter wavelengths. The pulsar has not been found 
at any radio wavelength. 

We found strong linear polarization in the di- 
rection towards the SNR. The electric field vec- 
tors are superimposed on continuum contours in 
Fig. 2. The mean fractional polarization is <~10%, 
calculated by dividing the polarized intensity im- 
age of the source by the total intensity image (both 
made with field no. 15 only). The total inten- 
sity values were truncated at the 3-ct level which 
is ^5mJybeam _1 before division, but all (bias 
corrected) polarized intensities within the source 
boundaries given by the total intensity cutoff were 
used. In the vector display, the vectors were trun- 
cated at the 3-cr level. We note that a polarization 
map using seven mosaic positions looked similar 
but was noisier because of the few percent radial 
instrumental polarization toward N157B from the 
surrounding fields. 

The electric vector alignment is fairly uniform 
but the polarized intensity is surprisingly patchy. 
The highest fractional polarization is 15.4% to- 
wards 05 h 37 m 44?5 and -69°10'07'/5 3 , whereas 
at the position of the brightest total intensity 
05 h 37 m 41 s 5 and _69°09'31?0 it is 5%, which is the 
lowest fractional polarization (excluding blanked 
regions). At the position of the X-ray peak at 
05 ?t 37 m 47?5 and -69°10'17'/5 the polarized inten- 
sity is 7%. 

3.2. 6 cm image 

The total intensity image of N157B at 6 cm 
(shown in Fig. 3) was made using eight available 
configurations of the Compact Array (see Table 1) 
and also has the synthesized beam of HPBW ~2". 
The fractional polarization towards N157B (calcu- 
lated using only field no. 15) is only ~3%, which 
is close to the noise level. We therefore detect 
no reliable polarization at 6 cm and conclude that 
there must be strong Faraday depolarization be- 
tween 3.5 and 6 cm. We note that at 6 cm the 
distribution of the closest antenna spacings repre- 
sents lower spatial frequencies than those seen at 
3.5 cm. If the polarization is all on a fine spatial 

3 all position are in J2000 




5 h 37 ra 56 B 52" 48 s 44 s 40 s 36' 

RA (J2000) 



Fig. 1. — Radio continuum grayscale and contour 
image of N157B at 3.5 cm with HPBW of 2'/67 x 
l'/80. The contour levels are: 1.1, 1.6, 3.3, 4.9, 6.5, 
8.2 and g.SmJybeam^ 1 . 




5 1 37 m 56 a 52 = 4B a 44 s 40 = 36' 

RA (J2000) 



Fig. 2. — Image of the total intensity of N157B 
at 3.5 cm with E-vectors superimposed. A vector 
length of 1" represents a linearly polarized inten- 
sity of 1.7mJybeam _1 . The contour levels are: 
1.1, 1.6, 3.3, 4.9, 6.5, 8.2 and g.SmJybeam- 1 . 
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scale we may find a lower value of the polariza- 
tion at 6 cm. Inspection of Fig. 2, however, shows 
that while there is fine structure a significant part 
of the 3.5-cm polarization is uniform on scales of 
over 10" and so should appear at 6 cm. Depolar- 
ization effects are obviously present. 

4. Discussion 

N157B is very close to 30 Doradus and it is 
also embedded in an H II region around the OB 
association LH99 (Lucke & Hodge 1970). Chu et 
al. (1992) concluded that the southern part of the 
SNR is obscured and possibly interacts with the 
dark cloud within the H II region. With our high 
radio sensitivity we can see that the SNR does ex- 
tend well to the south of the bright northern com- 
ponent. With a size of 30 pcx 18 pc, it is the largest 
Crab-type SNR known. The feature at 05 h 'S7 m 52 s 
and — 69°10'22" and the elongated structure cen- 
tered at about 05' l 37 m 50 s and -69°11'10" might 
indicate the beginning of a transition into a com- 
posite remnant like SNR 0540-693. The dense in- 
terstellar environment of the 30 Doradus region 
may be conducive to the formation of a shell and 
the irregular structure of N157B may be accounted 
for by a complex density distribution. 

Chu et al. (1992) found that N157B contributes 
up to 80% of the flux density in this region at 0.843 
GHz , the rest coming from the surrounding H II 
region. They suggested that the flat radio spectral 
index (a) of —0.1 (S v ~ v a ) indicated from earlier 
measurements (Milne et al. 1980) should be char- 
acteristic of the SNR itself. However, they also 
found that the optical emission lines from N157B 
originate in conditions associated with blast waves 
and they therefore concluded that the correspond- 
ing non-thermal radio emission should not have 
the flat spectrum typically associated with pulsar 
stimulation. Our results below differ somewhat 
from the previously published values derived from 
low resolution data but still show that N157B has 
a flat radio spectrum. 

Data used for determining the radio spectrum 
of N157B are listed in Table 2. To obtain a con- 
sistent background and source extent, the images 
were all convolved to the 43" x 45" resolution of 
the MOST image and the integrated flux density 
calculated over the same area in each convolved 
map. A background intensity around the source 



was determined from the average of several areas 
around the SNR but far enough from 30 Doradus 
to avoid contamination. 

The resultant integrated flux densities listed in 
Table 2 give a mean spectral index for N157B 
of —0.19 (see Fig. 4). Brightness-brightness plots 
also show that the spectrum does not vary signif- 
icantly across the whole extent of the SNR. This 
argues that the entire SNR is pulsar-powered. 

In the X-ray band, however, the spectrum of 
N157B differs slightly relative to other such rem- 
nants. It has a power-law energy slope of —1.5 
which is steeper than those of other Crab-like rem- 
nants with values of about —1.0 (Wang & Got- 
thclf 1998). The X-ray and entire spectrum of 
N157B also strikingly resemble that of SNR 0540- 
693, leaving no doubt of the plcronic origin of 
N157B. The SNR 0540-693 contains a pulsar and 
is a composite remnant, sharing the characteris- 
tics of pure shell remnants and pure plerions (Se- 
ward 1989; Manchester et al. 1993). This latter 
composite SNR is often called a Crab-type rem- 
nant, because the core is Crab-like and the faint 
shell component was only investigated later. This 
further similarity between N157B and 0540-693 is 
more evidence that N157B is beginning the tran- 
sition to a composite remnant. 

The polarization properties of the 30 Doradus 
region were previously investigated by Dickel et al. 
(1994) at longer wavelengths. They used ATCA 
observations at 13 cm but concluded that no polar- 
ized emission was detected in any region of N157B 
above the level of instrumental polarization (less 
than 2%). We report a significant level of polar- 
ized emission at 3.5 cm from N157B for the first 
time. This relatively high level of 10% is typical 
for a Crab-type remnant. The data at 6 cm give 
a marginal detection which we consider an upper 
limit of 3% fractional polarization at the peak of 
N157B. This low level of polarized emission at this 
wavelength implies very strong internal depolar- 
ization. The observed polarization could also be 
reduced by the thermal emission from the associ- 
ated H II region around LH99 which may be inter- 
twined with the SNR so that its thermal electrons 
may cause some of the Faraday depolarization as 
well. 

To estimate the internal rotation necessary 
to cause the observed depolarization, we use a 
slightly modified version of the relation given by 
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Table 2 
Flux Densities of N157B. 



Frequency (GHz) 


Flux Density (Jy) 


References 


0.843 


2.86 ±0.24 


Anne Green (1997, private comunication) 


1.400 


2.64 ±0.16 


Data of Mebold et al. (1997) 


2.378 


2.36 ±0.14 


Data of Dickel et al. (1994) 


4.740 


2.20 ±0.10 


Data of this paper 


8.600 


1.82 ±0.12 


Data of this paper 



Burn (1966): 

P(X 2 )/P(Xi)=sin(6)/6 (1) 

where P is the fractional polarization at A2, the 
longer wavelength, and Ai, the shorter wavelength; 
6 is the rotation in position angle of the measured 
electric vector between the two wavelengths. We 
will adopt the possible detection of 3% polariza- 
tion at 6 cm and so use a rough mean value for 
-P(A2)/-P(Ai) of 0.3 which requires an average ro- 
tation of 135°. Then using the standard formula 
for Faraday rotation we can estimate the line-of- 
sight magnetic field strength: 

5 = 7.9 x 10 5 \ 2 J N e Bxdl (2) 

where 5 is in radians, A is in meters, N e is the 
electron density in cm~ 3 , B is the magnetic field 
strength in Gauss, and / is the path length through 
the SNR in pc. Because the Faraday rotation 
increases progressively along the path, we must 
actually find the adopted rotation along half the 
depth along the line of sight. We do not, of course, 
know this depth directly but choose its value to be 
approximately equal to the short observed axis of 
the SNR or a half depth of 10 pc. From an anal- 
ysis of ROSAT and ASCA X-ray spectra, Wang 
& Gotthelf (1998) derived an electron density of 
0.6/~ 5 cm~ 3 , where / is the filling factor of the 
object. Because Crab-type SNRs should be rea- 
sonably uniform inside, we adopt a value for / of 
about 0.5 to get a mean density of about 2cm -3 . 
These numbers give a line-of-sight magnetic field 
strength of 30^iGauss. Most Crab-type SNRs 
have rather uniform magnetic fields; over the lim- 
ited region where it can be observed, N157B shows 



more variation than most. It may be that the 
field is directed mainly toward us, which would re- 
duce the observed polarization at all wavelengths 
and also give a large Faraday rotation. No cor- 
rection has been made for an unknown field ori- 
entation which would tend to increase the field 
strength in the remnant. We should emphasize 
that this result is only an order of magnitude es- 
timate. Neglected factors include the uncertain- 
ties and probably variations in the depolarization, 
magnetic field reversals, density uncertainties and 
variations. For example, significant clumping can 
change the apparent integrated electron density 
because the X-ray emission depends upon the den- 
sity squared but the Faraday rotation depends lin- 
early on the density and the sum of squares does 
not equal the square of sums. 

If we know the magnetic field strength, we can 
then determine the energy in relativistic electrons 
from the intensity of the synchrotron radiation. 
Ginzburg & Syrovatskii (1965) give the approxi- 
mate formula: 

E e = Y, v x R 2 x Axareaxpathlengthx B~ z/2 (3) 

where E e is the total energy in relativistic elec- 
trons, X„ is the surface brightness of the rem- 
nant at a given observing frequncy v, R is its 
distance, and A is a parameter that depends on 
the frequency, the spectral index and the adopted 
cutoff frequencies for the synchrotron radiation. 
There appears to be a break in the spectrum 
of the non-thermal emission between the radio 
and X-ray (Wang & Gotthelf 1998) but we do 
not know specifically where it occurs. We shall 
choose an upper frequency of 10 12 Hz and a lower 
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5 h 37 m 56 s 52 s 48 s 44 s 40 s 36 ! 

RA (J2000) 



Fig. 3. — Radio continuum grayscale and contour 
image of N157B at 6 cm with HPBW of 1"79 x 
l'/74. The contour levels are: 0.9, 1.3, 2.4, 4.0, 
5.4, 8.1, 10.8 and lS.SmJybeam- 1 . 



io-o 



5.0 - 




1.0 - 

0.5 I I I I I I i i ' 

0.5 1.0 2.0 5.0 10.0 

FREQUENCY (GHz) 

Fig. 4. — Radio spectrum of N157B. The slope of 
the fitted line is —0.19. 



one of 10 7 Hz, with a spectral index of —0.19, 
A = 3.7 x 10 19 at 8.6 GHz where the flux density is 
1.82 Jy. With the sizes and magnetic field strength 
given above, we find that the mean relativistic 
electron energy density is 4 x 10 -10 ergscm~ 3 . 
This is about 10 times the magnetic energy den- 
sity. While within the errors, we cannot definitely 
state that the values are different. It is interest- 
ing to note that in the three objects for which we 
have now measured the relative energies by this 
method, each has between 3 and 10 times greater 
relativistic electron energy than magnetic energy 
(Kepler's SNR Matsui et al. (1984) and N23 SNR 
Dickel & Milne (1998)). 

Wang & Gotthelf (1998) give an approximate 
temperature for the thermal gas in N157B of 0.4 
- 0.7 keV. Choosing 0.6 keV and the approximate 
density of 2 cm -3 used above, the thermal energy 
density is about 3 x 10 -9 ergcm~ 3 . Thus the ther- 
mal energy is about 10 times that in relativistic 
electrons and 100 times that in the magnetic field. 
Even in this Crab-type SNR the thermal energy 
dominates. 

In summary, the filled structure and uniform 
spectral index of N157B are Crab-like but it also 
has some characteristics similar to those of the 
composite SNR 0540-693. Although the dense, 
irregular surroundings may be partly responsible 
for the non-uniform brightness and outline of this 
SNR particularly to the south, the overall mor- 
phology and spectral index suggest that this SNR 
is just beginning the transition from a pure Crab- 
like remnant to a composite remnant. It will be 
interesting to see how it changes at radio, optical 
and X-ray wavelengths in the future. 

We have detected for the first time a signifi- 
cant level of linearly polarized emission in this ob- 
ject at 3.5 cm wavelength. The encompassing H 
II region and clumpy internal structure within the 
SNR may cause sufficient Faraday rotation to re- 
duce the polarization to below detectable levels 
at 6 cm and make the large variations in the po- 
larization observed at 3.5 cm at different positions 
around the remnant. Observations at higher fre- 
quencies, e.g. 23 GHz, are needed to enable more 
accurate determination of the Faraday rotation. 

We thank M. Costa, D. Milne and H. Dickel 
for participating in the observing. We thank L. 
Staveley-Smith for the ATCA image at 21 cm and 
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